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On the Mott glass in the one-dimensional half-filled charge density waves 
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We study the effect of impurity pinning on a one-dimensional half-filled electron system, which 
is expressed in terms of a phase Hamiltonian with the charge degree of freedom. Within the 
classical treatment, the pinned state is examined numerically. The Mott glass, which has been 
pointed out by Orignac et al. [Phys. Rev. Lett 83 (1999)2378], appears in the intermediate 
region where the impurity potential competes with the commensurate potential. Such a state is 
verified by calculating the soliton formation energy, the local restoring force around the pinned 
state and the optical conductivity. 
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§1. Introduction 

In one-dimensional electron systems, Anderson localization corresponding to the impurity pinning 
takes place even for weak impurities. 1-3 ) The pinning is enhanced for repulsive interactions, as 
found in the charge density waves 4, 5 ) where the spatial variation of the phase spreads widely due 
to randomness. 6 ^ However, such a pinning becomes complicated for the system with both the 
commensurate band filling and the repulsive interaction, in which the phase is also pinned by the 
commensurate potential, 7 ) i.e, umklapp scattering. It has been maintained that the Mott insulator 
with a charge excitation gap occurs for large commensurate potential and the compressive Anderson 
glass with no excitation gap appears for the weak random impurities. 8, 9 ) Such a competition has 
been further investigated by Orignac, Giamarchi and Doussal (OGD), 10 ) who applied the variational 
method to the replicas for the impurity potential. For the case of the impurity potential being 
much larger than the commensurate potential, one obtains the Anderson glass, in which the optical 
conductivity exhibits (t(uj) oc lu 2 i.e., with no significant effect of the commensurate potential. 
The opposite case leads to the Mott insulator where the compressibility becomes zero and the 
gap appears in the optical conductivity. However there exists the intermediate case called the 
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Mott glass where the compressibility becomes zero but the conductivity shows the gapless behavior 
similar to that of the Anderson glass. 10, n ) Such a state comes from two kinds of length scales ; one 
is the localization length and the other is a width of the kink followed by the soliton formation. 
The vanishing of the compressibility indicates a finite energy to create the soliton of the phase 
of the charge density wave, i.e., the relevance of the commensurate potential. Thus the energy 
for the soliton is determined mainly by the local region around the kink. On the other hand, the 
absence of the gap in the conductivity comes from the property of the impurity pinning with the 
localization length, which is different from that of the kink. Thus these two kinds of pinning could 
be compatible and lead to a certain region of coexistence. 

Such a coexistence is explored in the present paper by examining the properties of the ground state 
and the excited state of a system, which lead to the competition between the impurity potential 
and the commensurate potential. We study numerically a half-filled electron system with repulsive 
interaction based on the phase Hamiltonian. For the simplicity, we focus on the charge degree of 
freedom. Based on our preliminary work, 12 ) we demonstrate that such a Mott glass state does 
appear with increasing the commensurate potential as the crossover from the impurity pinning to 
the commensurability pinning. In §2, formulation is given and the pinned states are examined to 
determine the competition between the Anderson glass and the Mott glass. In §3, the competition 
between the Mott glass and the Mott insulator is examined by calculating the local fluctuation 
around the pinned state, e.g., the restoring force and the optical conductivity. Discussion is given 
in §4. 

§2. Model and the Pinned State 

The formulation for calculating a model numerically is given within the classical treatment. 
The parameter for the competition between the impurity potential and commensurate potential is 
derived. The excitation energy for the soliton formation is calculated in the presence of impurity 
where the appearance of the finite energy is identified with a crossover from the Anderson glass to 
the Mott glass. 

2.1 Model 

We consider a one-dimensional Hamiltonian with a length L given by 4 ) 

N imp 



H = I dx 



o 



A i^—J - V cos ( 9 + 2k F x)8(x - Rj) - B cos(2fl) 



(1) 



where 9(x) denotes a phase variable of the charge density wave. 13 ) The V term denotes the impurity 
potential where Ni mp is the total number of the impurity and Rj is the location of the impurity at 
the j-th site. Equation (1) without the V term denotes the Hubbard model with a half-filled band 
and repulsive interaction, C/, 14 ) in which the quantum fluctuation and the spin degree of freedom are 
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discarded. Equation (1) describes the classical charge density wave where A = vf/(4tt), B ~ U/a 
with Fermi velocity, vp, and the lattice constant, a, 13 ) The A term leads to the uniform 6(x) while 
the V term and the B term lead to spatially varying 6(x). We consider the A term is much 
larger than the V term and the B term, i.e., ANi mp /L » V and A(Ni mp /L) 2 » B. In such a 
case, the length corresponding to the impurity concentration is much smaller compared with the 
characteristic length determined by V and B terms. By taking the average spacing of the impurity, 
L/N imp , as a new lattice constant, eq. (2) is rewritten as the following discretized model. 15 ) 



i=i L 



B 

A ni (9 j+l - 8j) 2 - Vcos(9j - Q) cos(2%) 



(2) 



where rii = Ni mp /L and Cj( = ~ 2/cf-Rj) is a random variable with a Fermi momentum kp. The 
characteristic length Iq, which denotes the pinning by the impurity potential, is determined by 
the following variational method. Considering eq.(2) with B = 0, one obtains the energy gain of 
the second term given by —V/y/njlo per lattice site, but the increase of the first term given by 
~ A/nily. Minimizing the total energy consisting of these two quantities, one obtains Iq as 6 ) 

where Kq = 7r 2 /3 and 

V 



" AirAm ' ^ 
The present paper treats the case of small e(<< 1) corresponding to the weak impurity pinning. 
On the other hand, eq. (2) with V = gives the excitation followed by the soliton formation with 
an energy, 4v / 2(AB) 1 / 2 , and a characteristic length d, 

d= 4) ■ (5) 



B J 

Since Iq is the pinning length for the V term and d is the width of the soliton for the B term, the 
competition between these two terms is characterized by the quantity 7 ) 

k= l ° i • (6) 
d(2K )i 

It is expected that the crossover from the impurity pinning to the commensurability pinning occurs 
for l ~ o(l). 

For the convenience of the numerical calculation, eq. (2) is rewritten as 

H = AriiJ2 [(°3+i ~ e if ~ 4vrecos(% - Q - /3cos(2^-)] , (7) 
i=i 

where 

B 2 

An i K§ 
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Equation (7) shows that the energy is measured in the unit of Ani and the lattice constant is taken 
as 1/nj. By minimizing eq. (7) with respect to Oj, we obtain the equation for Oj as 



where j = 1, ■ • • , N. A periodic boundary condition is taken where Oq = and 0^ imv +i = 0±. 

The following calculation is performed by taking two kinds of parameters, e and Iq which corre- 
spond to the variation of V and B respectively. 

2.2 spatial variation of Oj 

We calculate eq. (9) iteratively with several choices of initial values and stop the calculation 
when the absolute value of the l.h.s. becomes less than 10 -7 . The strength of the impurity pinning 
is chosen as e = 0.01, which denotes the weak pinning and leads to the pinning length, u^Iq ~ 22. 
In this case, considering the periodic boundary condition, one expects that there are about ten 
regions of Oj exhibiting the rapid variation of the order of ir for a system with Ni mp = 400. 

In Fig. 1, we show the numerical results for Ni mp = 400, e = 0.01, Iq = 0, 0.5 and 1.5 respectively. 
There are many solutions corresponding to the metastable states, which have the energy close to 
the ground state. Figure 1 may be regarded as the ground state since the corresponding energy is 
the lowest one within the present calculation. The variation of 0j becomes small with increasing 
Iq i.e., the commensurability energy. Such a state is the solution for the uniform sector since the 
spatial variation becomes constant (e.g. 0j = 0) in the limit of large Iq. 

Figure 2 shows the spatial variation of 0j for the excited state in which the soliton exists. For 
one of the the kinks is found at Xj ~ 170 although the initial value of the kink is taken 
at Xj = 200. This indicates that the kink is determined in order to gain the energy from the 
impurity potential even for the large Iq. The kink begins to move from Xj ~ 170 to Xj ~ 140 
with decreasing Iq while the other kink also moves at lower Iq to gain the energy from the impurity 
potential. However the location is fixed for the case of Iq ^ 0.6 indicating a crossover from the 
commensurability pinning to the impurity pinning. Such a state is regarded as the solution from 
the soliton sector since the spatial variation reduces to the solution of the soliton-antisoliton in the 
limit of large Iq. 

2.3 excitation energy 

There are two kinds of excitations from the ground state. One of them is the excitation within 
the uniform sector, which is gapless one due to the continuous distribution of the energy. The 
other is the excitation followed by the soliton formation where the lowest one is a state with the 
soliton-antisoliton. The energy difference between these two states, 5E exc , is calculated by 



{26 j - j+1 - 6j-i) + 4vresin(^ - Q + 2/3sin(2^) = , 



(9) 



$E exc — SE so i o~Eg roun( [ . 



(10) 



4 







200 



400 



Fig. 1. Spatial variation of 0j for Ni mp = 400, e = 0.01 with fixed lo =0 (solid line), 0.5 (dashed line) and 1.5 
(dotted line). The unit of Xj(= 1, • • • Nim P ) is given by l/n». 



Quantities 5E so i and 6E groun d, which denote the deviation of the the energy from that of the state 
with 9j = 0, are calculated by substituting 9j of the respective sector into the following equation: 



where E denotes an expectation value of H given by eq. (1). 

It should be noted that there is a critical value ls(— 0.6) as shown by the arrow in Fig. 3 where 
5E exc (> 0) is very small and almost constant for Iq < 1$. Although the spatial variation of the 
soliton sector is apparently retained even for small Iq, the effect of the commensurability energy is 
negligibly small and then such a spatial variation can be regarded as one of the metastable states 
in the absence of the commensurability energy. For Iq ^ Is, one obtains SE exc ~ 0.25/o — 0.8, which 
indicates that the formation energy of soliton-antisoliton is finite but is reduced by a finite energy 
due to the energy gain by the impurity potential. 

§3. Fluctuation around the Pinned State 

In order to find the crossover from the Mott glass to the Mott insulator, the degree of the local 
charge fluctuation around the ground state is examined by calculating the response to the external 
field, which gives rise to the optical conductivity. 



5E = (- + B)/(3AKq/1 2 ) 




(11) 
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Fig. 2. Spatial variation of 9j for Ni mp = 400, e = 0.01 for the excited state followed by a formation of a soliton- 
antisoliton where lo ~0 (solid line), 0.5 (dashed line) and 1.5 (dotted line). 



3. 1 local pinning frequency Uj 

We examine the restoring force for 8j at the j-th site when the small deviation of 6j occurs from 
the equilibrium value. The magnitude of the force depends on the location of the lattice site since 
the pinning comes from not only the commensurate potential but also the impurity potential. Such 
a quantity is estimated from the increase of the energy by adding the external electric field E' 
which gives an additional term, — E0, to the r.h.s. of eq.(l) where E = eE f /ir with e being the 
electric charge. 16 ) In this case, the lattice model is rewritten as 



N imp 



H = An, t J2 [(0. 



'i+i 



Ojf - 47recos(0j - (j) - /3cos(2^) - jE9j 



(12) 



3=1 



where E = l 2 /(3AK ) and 



7 = 



3(vre)3 



(13) 




The resultant equation for 0j is given by 



(20 j - j+ i - ^_i) + aam(6j - Q) + 2/3sin(2^) - ~/E = . 



(14) 



By using 9j(E) of the solution of eq. (14), we calculate 



in 



hm (e 3 {E)-90))hE . 



(15) 
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Fig. 3. Excited energy SE exc as a function of lo for Ni mv = 400 and e = 0.01, which is averaged over 400 samples. 
The arrow Is denotes a critical value above which the soliton formation is followed by a finite energy. 



The quantity uj is the local pinning frequency since the increase of the energy in the presence of 
the small external field can be expressed as Uj(6j(E) — 6*j(0)) 2 /2 , which is calculated by expanding 
eq. (12) in terms of 9j(E) - 9(0). 

In Fig. 4, lo dependence of < Uj > is shown for Ni mp = 4000 and e = 0.01 where <> denotes 
the spatial average and the closed square is the result averaged over 20 samples. The calculation 
of ojj has been performed by choosing E = 0.001, which is small enough to calculate eq. (15). 
For Iq < 0.4, the quantity < Uj > is constant indicating no visible effect of the commensurability 
pinning. The commensurability pinning dominates for Zo~0.6, in which Uj oc Iq. It seems that a 
crossover from the impurity pinning to the commensurability pinning occurs in the interval region 
of 0.4<I^<0.6. 

3.2 root mean square of u>j 

The property of the spatial variation of ujj is examined by calculating the root mean square of 
Uj . The normalized root mean square of ujj is given by 

z = ±_J({ UJ] -{u J] )Y) . (16) 

{Uj) 

The larger z indicates the broader distribution of the pinning force. It is expected that z of the 
impurity pinning is much larger than that of the commensurability pinning. The fluctuation of z is 
large for the choice of the present sample with Ni mp where z s is used to distinguish the respective 
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Fig. 4. lo dependence of < ujj > for Ni mp — 4000 and e = 0.01, which is averaged over 20 samples. The limiting 
behavior is evaluated by taking E — 0.001. The error bar becomes negligibly small in the visible scale with 
increasing lo- 



sample. Thus, we further examine the root mean square of z s given by 



Figure 5 (dashed line and closed circle ) shows the normalized root mean square of < Uj >, i.e. 

< z s > s for Ni m p = 400, e = 0.01 and E = 0.001 where 400 samples are used for s. For Iq % 0.6, 

< z s > s is almost constant while it decreases for Iq ^ 0.6 due to the increase of the commensurability 
pinning. The result for the soliton sector is shown by the open circle where it remains finite even 
for large Iq. This is understood as follows. When the impurity is absent, the location of the kink 
is easy to move under the external field resulting in the extremely small ojj . However the presence 
of the impurity pinning prevents such a movement leading to finite Uj, which gives the variation 
similar to the impurity pinning. The deviation from that of the uniform sector is found for Iq ^ 0.7, 
which is slightly larger than 1$. This indicates a possibility that the fluctuation around the pinned 
state even for Iq ^ Is still resembles that of Iq = 0. The inset shows the root mean square of such 
z s which gives a critical value Iq ~ 1.0, above which the spatial fluctuation of Uj decreases rapidly 
indicating a crossover into the Mott insulating state. 




(17) 
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Fig. 5. lo dependence of the normalized root mean square of uij for Ni mp — 400, e = 0.01 and E — 0.001 which is 
averaged over 400 samples. The closed circle and solid line (open circle and dotted line) denote the results for the 
ground state (the excited state with soliton) . The inset denotes the corresponding root mean square of z s given by 



This tendency can be also seen from the distribution of ujj , which is defined by 

PH = T r - E ^UJ-Uj) . (18) 

In Fig. 6, lo(= ojj) dependence of P is shown for JVj mp = 4000, e = 0.01 and E = 0.001 where Zo = 
(closed circle), 0.5 (open circle), 0.8 (closed triangle) and 1.5 (closed diamond). The average over 
20 samples are used for the numerical accuracy of small uj. The quantity P with Iq = is almost 
the same as P with Zo = 0.5 and still resembles P with Zo = 0.8 while it is much different from P 
with Zo = 1.5. This indicates that the distribution changes the property for Zo > Ig > where Iq is 
larger than at least 0.8. 



3.3 optical conductivity 

Based on Figs. 3, 4, 5 and 6, it is reasonable to admit that the intermediate region with Is <h < 
Iq may correspond to that of the Mott glass. However the estimation for l& is rather complicated 
compared with lg- Then we examine the optical conductivity since the Mott glass is maintained to 
exhibit the optical conductivity given by a (to) oc aA 10 ) We calculate the conductivity within the 
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0)/<0)j> 



Fig. 6. Normalized distribution of ujj, P{u)/ < uJj >), for Ni mp — 4000, e = 0.01 and E = 0.001 where lo — (closed 
circle), 0.5 (open circle), 0.8 (closed triangle) and 1.5 (closed diamond). The results are obtained by averaging over 
20 samples. 



classical treatment, by noting the fact that the electric dipole moment at the j'-th site is given by 
(2kF/ir)eyj = (e/ir)9j with yj being the displacement. 16 ' 1 In the presence of the electric field Eoe lu}t , 
the equation of motion for the charge density wave is given by 

*d 2 yj m* dyj „ i^t * 2 /m\ 
m -pt = — — -± + eE e - m , (19) 

where m* is the effective mass and the life time, r, is introduced for the convergence of the numerical 
calculation. The quantity loj, which corresponds to the restoring force, m*uJjyj is calculated from 
eq. (15). The quantity yj is calculated as 

eEo/(iuim*) 

Vi = 7 ~J\ > ( 20 ) 



which leads to the local conductivity at the j-th site 



Thus one obtains the real part of the total conductivity as 



2k F e 2 r . 

(Jj{uj) = -, ^- . (21) 



n "imp -. 

a ^ = ^v— E 7 — —1 ' (22) 
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where the conductivity is normalized such that the summation of ct(uj) with respect to u becomes 
equal to unity. In Fig. 7, optical conductivity is shown for Ni mp = 4000, e = 0.01 and r = 20 
where respective results are for Iq = (closed circle), 0.5 (open circle), 0.8 (closed ) and 1.5 (closed 
diamond). The conductivity for Iq = denotes the conventional impurity pinning of the charge 




CO 



Fig. 7. Normalized optical conductivity for N imp — 400, e — 0.01 and r = 20, which are averaged over 100 samples. 
Respective results are for lo = (closed circle), 0.5 (open circle), 0.8 (closed ) and 1.5 (closed diamond). The 
dashed curve corresponds to the conductivity of the Drude's formula (i.e., ujj — » 0). 



density wave where the frequency of uj corresponding to the peak is essentially the same as < Uj >, 
i.e., the average of Uj shown in Fig. 4. Such a frequency increases with increasing lo, which is 
consistent with the result of Fig. 4. The width of uj is determined by the root mean square of 
Uj, (Fig. 5). The width for small lo, which is determined by the distribution of the pinning force, 
decreases for large Iq due to the increase of the commensurability pinning. It is found that the 
peak of uj is slightly reduced in region of the Mott glass (Is < Iq < Iq), since the width of <t(u) 
in the Mott glass becomes large due to the competition between the impurity pinning and the 
commensurability pinning. For comparison, the conductivity of the Drude's formula is shown by 
the dashed curve where the width is chosen as r = 20. The conductivity for Iq =0, 0.5 and 0.8 
exhibits <r(u) oc uj 2 for small uj within the numerical accuracy of the present calculation indicating 
the property of the Anderson glass. In order to examine such a region, we calculate a coefficient, 
A, of a(uj) which is given in the limit of small u by 

ct(uj) ~ Auj 2 . (23) 
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In Fig. 8, the quantity A is shown by taking oj — > 0.01, where the dashed line denotes the guide for 
the eye. The inset depicts the procedure to calculate A by an extrapolation of A(r) for 1/r = 1, 
0.1 and 0.05 since eq. (23) depends on r. With increasing lo, the quantity A stays almost constatnt 
for small Iq while it decreases rapidly for Iq ~ 1.0. However, the data for Iq > 1.0 is beyond the 
present calculation, which is based on the model of the classical charge density wave, and eq. (22) 
showing the gradual decrease of A even in the Mott insulating phase. Thus we concludes Iq — 1-0, 
which is shown by the arrow in the main figure. 




Fig. 8. lo dependence of A, which denotes <t(oj) ~ Aui 2 in the limit of small lo. In the inset, A is obtained from the 
value of the vertical axis where the extrapolation for 1/r is shown for lo = 0.8. 1.0 and 1.2. The steep decrease of 
A leads to Ig — 1-0 as shown by the arrow. 



§4. Discussion 

We have examined the competition between the impurity pining and the commensurability pin- 
ning of the charge density wave, which leads to three kinds of phases as a crossover. It is found that 
the Mott glass state exists in the intermediate region of 1$ < lo < Ig with Is — 0.6 and Iq — 1-0 
while there are the region of the Anderson glass for lo < Is and that of the Mott insulator for 
h <To- 

The intermediate region corresponding to the Mott glass seems to be consistent with that of 
OGD, 10 ) who obtained 0.98 < d/lo < 1.58 in their notations. The Mott glass undergoes the 
crossover to the Anderson glass (or Mott insulator) with decreasing ( or increasing) lo due to the 
effect of the impurity, which has the wide distribution of the pinning force (i.e. restoring force). 
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Then there is a tiny domain of the soliton formation in the state of the Mott glass. Even in the 
presence of such a state, the optical conductivity similar to the Anderson glass could be expected 
due to the contribution from the slowly varying region 17 ) since the conductivity originates in the 
distribution of the restoring force in the uniform sector. The conductivity with low frequency 
resembles that of the Anderson glass, which denotes no gap in the conductivity, i.e., <j(oj) oc lo 2 
for small u. Note that it is beyond the present calculation to examine if the conductivity exhibits 
a property a(uj) oc o; 2 ln 2 (o;r). 1 ' 2 ) The existence of the gap for the formation of the phase soliton, 
which leads to the variation of the electron number, is equivalent to the absence of compressibility 
since a finite energy is needed to add or remove an electron. However the case of doping of electron 
or hole, which corresponds to the ground state with the soliton formation, leads to the absence of 
the Mott insulator and the optical conductivity with no gap as expected from the soliton sector of 
Fig. 5. 

Here we comment on the present mode of the charge density wave, which is relevant to the 
electron system. The latter case includes two additional factors, i.e, the quantum fluctuation and 
the spin fluctuation. These factors have at least the effect of reducing the magnitude of V and B 
in eq. (2), which could be taken into account by adjusting parameters. However the novel effect 
can be expect when the depinning begins. This problem will be discussed in a separate paper. 
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